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The newly created Bose-Einstein condensates (BEC) of weakly interacting alkali-metal 
atoms [1] stimulated a large number of theoretical investigations (see recent reviews [2]). 
Most of these works are based on the assumption that the properties of BEC are well de- 
scribed by the Gross-Pitaevskii (GP) mean-field theory [3]. The validity of the GP equation 
is nearly universally accepted. 

The experimental realization of quasi-one dimensional (ID) and quasi-two dimensional 
(2D) trapped gases [4-6] stimulated many theoretical interests. The theoretical aspects of 
the BEC in quasi-lD and quasi-2D traps have been reported in many papers [7-17]. For 
the case of dimensions d < 3, it is known that the quantum-mechanical two-body t-matrix 
vanishes [18] at low energies. Therefore the replacement of the two-body interaction by the 
t-matrix, as done in deriving the GP mean-field theory, is not correct in general for d < 3 



The density functional theory (DFT), originally developed for interacting systems of 
fermions [20], provides a rigorous alternative approach to the interacting inhomogeneous 
Bose gases. The main goal of this letter is to develop a Kohn-Sham (KS) like time-dependent 
theory for bosons. 

We consider a system of N interacting bosons in a trap potential V ext . Assuming that 
our system is in a local thermal equilibrium at each position f with the local energy per 
particle e(n) (e is the ground state energy per particle of the homogeneous system and n is 
the density), we can write a zero temperature classical hydrodynamics equation as [8] 



[12,19]. 
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where v is the velocity field. 

Adding the kinetic energy pressure term, we have 
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We define the density of the system as 



n(f,t) =| #(f,t) 



(4) 



and the velocity field v as 
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From Eqs.(l) and (5), we obtain the following KS like time-dependent equation 

h 2 „ 2 T t, t 9(ne(n)) T 

If the trap potential, V ex t, is independent of time, one can write the ground-state wave func- 
tion as \I/(r, t) = $(r) exp(— i/xt/h), where fx is the chemical potential, and $ is normalized 
to the total number of particles, J df | $ | 2 = N. Then Eq.(6) becomes [21] 

where the solution of the equation (7) minimizes the KS energy functional in the local density 
approximation (LDA) 

h 2 

E = N<$\ -—V 2 + V ext + e(n) | $ >, (8) 
2m 

and the chemical potential fx is given by fx — dE/dN. Eq.(7) has the form of the KS 
equation. 

The ground-state energy per particle of the homogeneous system e(n) for dilute 3D [22] 
and dilute 2D [23] Bose gases are 

e(n) = ^-a 3D n[l + ^=(nal D fl 2 + 8(y - V3)na 3 3D Hna 3 3D ) + ...], (9) 

and 

e(n) = | Hna 2 2D ) ^ (1 + | \n(na 2 2D ) \' 1/5 ), (10) 

m 

where a 3 D and a<iD are 3D and 2D scattering lengths respectively. 

For ID Bose gas interacting via a repulsive 5-function potential, g5(x), e(n) is given by 

[24] 

e(n) = ^ 2 e(7), (11) 

where 7 = mg/(h 2 n) and 

e(7) = ^3 y_ g{x)x dx, 

1 . „ q(x)dx . ,^ nS 

9 ( y) = -( 1 + 2A/_ i A 4' ( j_ g)2 ). (12) 
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For small values of 7, the following expression for e(n), 



<»,-!(.-£,/=*,+... m 

is adequate up to approximately 7 = 2 [24]. 
For a large coupling strength g [24] 

. h 2 n 2 n 2 2k 2 n._ 2 „ „. 

6m mg 

Eq.(14) is accurate to 1% for 7 > 10 [24]. 

In the limit of large N, by neglecting the kinetic energy term in the KS equation (7), we 
obtain an equation corresponding to the Thomas-Fermi (TF) approximation 

d(ne(n)) , . 

V ext + V d K n " = fi (15) 

in the region where n(r) is positive and n(r) = outside this region. 

Now we turn our attention to elementary excitations, corresponding to small oscillations 
of ty(r,t) around the ground state. Elementary excitations can be obtained by standard 
linear response analysis [25,26] of the equation (6), as resonances in the linear response. We 
add a weak sinusoidal perturbation to the time-dependent equation (6) 

= ( -2^ V2 + Vext + 9J ^ 1 + /+e ~" + f ~ e ^ (16) 
and assume that solution of Eq.(16) has the following form 

tf(f,t) = e -^ t/h mr) + u{r)e- iwt + v*(r)e iujt ], (17) 

where $(r) is the ground-state solution of Eq.(7). 

Linearization in small amplitudes u and v yields inhomogeneous equations 



on 2 



(18) 



where n—\ $(r) I 2 and 



h 2 „ 9 d(ne(n)) d 2 (ne(n)) , 
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Setting f± to zero in Eq.(18), we obtain coupled equations 

on 2 



d 2 (ne(n))^ 2 

Lv H V ; $ it = -fiwu, 

an" 1 



(20) 



which can be used to calculate the energies £ = Tiu of the elementary excitations. 

Now we describe the application of the time- dependent equation (6) to the case of non- 
linear dynamics. We turn to the limit of very strong coupling between the interacting bosons 
in ID, the so-called To nks- Girardeau gas [27]. In this impenetrable bosons case, the energy 
density e(n), Eq.(13), reduces to e(n) = h 2 Tc 2 n 2 /6m, and Eq.(6) reads [12] 

. , h 2 d 2 h 2 n 2 I 1 4\ t 

• n IF = '"^a? + v ~> + "1ST ' 1 >*• (21) 

with | tt(a;,t) | 2 dx = iV. 

For a general time- dependent harmonic trap V ext = mu 2 (t)x 2 /2, with the initial condition 
^/(x, 0) = $(a;), where $(a;) is the ground-state solution of the time-independent equation 

we show that Eq.(21) reduces to the ordinary differential equations which can provide the 
exact solution of Eq.(21). 

Indeed, if we assume that the solution, t), can be expressed as 

nx,t) = ^lMe-^^, (23) 

we obtain the following equations for A and (3 after inserting Eq.(23) into Eq.(21): 

A + u 2 {t)\ = A(0) = 1, A(0) = 0, (24) 

= ^, = °- ( 25 ) 
Thus, the ordinary differential equations Eqs.(22), (24), and (25) give the exact solution of 
Eq.(21), and the evolution of the density can be written exactly as 

n(x,t) = (1/X(t))n(x/X(t),0). (26) 
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For the case of free expansion, the confining potential is switched off at t = and the 
atoms fly away. In this case, Eqs.(24) and (25) can be integrated analytically leading to the 
following solutions for A and (3: 

X(t) = ^/l+cu 2 (0)t 2 , /3(t) = j-^-j arctan(cc>(0)t). (27) 

We note that self-similar solutions [28] of Eq.(21) were briefly discussed in Ref.[29]. 
In the large N limit, where the kinetic energy term in Eq.(22) is dropped altogether (the 
so-called Thomas- Fermi limit), the corresponding density is 

^= ~ wF 6 ^ ~ w?' (28) 

and for the Fourier transform n(k, t) = ^= f+™ n(x, t)e lkx dx we have 

n T F(k,t) = = Si v , 29 

where X(t) = [/i/(mc<j(0))] 1 / 2 A(t) and J\ is the Bessel function of the first order. 

The exact many-body wave function, ty B (xi,x 2 , ...,xw,t), of a system of N impenetrable 
bosons in a time dependent ID harmonic trap, can be found from the Fermi-Bose mapping 
[15] 

I ^b(xi,x 2 , ...x N ,t) 1 = 1 V F (xi,x 2 ,...x N ,t) I, (30) 
where typ is the fermionic solution of the time- dependent many-body Schrodinger equation 

= (31) 

dt 2m dx i 2 

with initial condition ^f(xi, x 2 , ...xn, 0) = $f(xi, x 2 , --Xn), where $f(xi, x 2 , ..xn) is the 
fermionic ground-state solution of the time-independent Schrodinger equation 

Therefore, for the exact density n B (x,t) = dx 2 ... dx N | * s ( x, x 2 , ...Xn, t) | 2 , we 
have 

1 N - 1 x 

n B (x,t) = J - T IHr-^l 2 , (32) 
A(tJ i=0 A{t) 

where 4>i(x) = Ci exp(— x 2 /2)Hi(x), Ci = 7r _1 / 4 (2 i i!) -1 / 2 , and Hi(x) are Hermite polynomials. 
Note that the evolution of n B (x,t) can be written as Eq.(26), corresponding to a time- 
dependent dilatation of length scale. 
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From the knowledge of n B (x,t) and riTF(x,t) one can evaluate the radii r(t) = 
n B (x, t)x 2 dxfl 2 and r TF (t) = n TF (a;, t)x 2 dxf/ 2 and ratio r(t)/r TF (t). This quan- 
tity is equal to 1 at any t for any iV. This circumstance explains why for a harmonic trap the 
ground-state density profile from Eq.(21) agrees well with the many-body results for systems 
with a rather small number of atoms iV m 10 [12]. As for a general trap potential we expect 
such agreement for much larger N. 

Using the following relation [30] 

J2 {2 m m)r l [H m {x)} 2 = (2 n+1 n\)- 1 {[H n+1 (x)} 2 - H n (x)H n+2 (x)}, (33) 

m=0 

we obtain an analytical formula for exact density t) 

n B (x,t) = ^c 2 N _ ie ^y^H[H N (x/X(t))] 2 ~ H N ^(x/X(t))H N+l (x/X(t))}- (34) 
Then the Fourier transform is given by 

MM) = '.-'^M^) 4- l( ^)], (35) 



where are Laguerre polynomials. Using the asymptotic formula of Hilb's type for the 
Laguerre polynomial [30], we have the asymptotic behavior of riB(k,t) as iV — > oo 

„ . N 2J 1 (V2N\(t)k) ^ /Arl/4x 

n ^-w, mxw + ° (Nl) ' (36) 

which is valid uniformly in any bounded region of k\(t). Eq.(36) for the case of t — 
is a rigorous justification of the Thomas- Fermi approximation [13,31] to a system of non- 
interacting ID spinless fermions in harmonic trapping potentials. 

Comparison of Eq.(36) with Eq.(29) shows that in the large iV limit the KS like time- 
dependent theory for ID impenetrable bosons in a time-dependent harmonic trap, Eq.(29), 
gives the same result as the exact many-body treatment, Eq.(36). Hence, we have rigor- 
ously proved that Eq.(29) correctly describes properties of ID Bose gas in a time-dependent 
harmonic trap in the limit of large N. 

In conclusion, we have developed a time- dependent KS like theory for bosons in three and 
lower-dimensional traps and have obtained coupled equations which can be used to calcu- 
late the energies of elementary excitations. For one- dimensional condensate of impenetrable 
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bosons in a general time-dependent harmonic trap, it is shown that the corresponding equa- 
tion reduces to the ordinary differential equations and gives the same results as the exact 
many-body treatment in the large N limit. 

We thank B. Tanatar for his interest and comments and E.B. Kolomeisky for informing 
us about Ref. [29] and for useful suggestions. 
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